The genomic landscape of acute lymphoblastic leukemia with intrachromosomal amplification of chromosome 21.

Intrachromosomal amplification of chromosome 21 defines a subtype of high-risk childhood acute lymphoblastic leukemia (iAMP21-ALL) characterized by copy number changes and complex rearrangements of chromosome 21. The genomic basis of iAMP21-ALL and the pathogenic role of the region of amplification of chromosome 21 to leukemogenesis remains incompletely understood. In this study, using integrated whole genome and transcriptome sequencing of 124 patients with iAMP21-ALL, including rare cases arising in the context of constitutional chromosomal aberrations, we identified subgroups of iAMP21-ALL based on the patterns of copy number alteration and structural variation. This large data set enabled formal delineation of a 7.8 Mb common region of amplification harboring 71 genes, 43 of which were differentially expressed compared with non-iAMP21-ALL ones, including multiple genes implicated in the pathogenesis of acute leukemia (CHAF1B, DYRK1A, ERG, HMGN1, and RUNX1). Using multimodal single-cell genomic profiling, including single-cell whole genome sequencing of 2 cases, we documented clonal heterogeneity and genomic evolution, demonstrating that the acquisition of the iAMP21 chromosome is an early event that may undergo progressive amplification during disease ontogeny. We show that UV-mutational signatures and high mutation load are characteristic secondary genetic features. Although the genomic alterations of chromosome 21 are variable, these integrated genomic analyses and demonstration of an extended common minimal region of amplification broaden the definition of iAMP21-ALL for more precise diagnosis using cytogenetic or genomic methods to inform clinical management.

Disruption to the FOXO-PRDM1 axis resulting from deletions of chromosome 6 in acute lymphoblastic leukaemia.

A common problem in the study of human malignancy is the elucidation of cancer driver mechanisms associated with recurrent deletion of regions containing multiple genes. Taking B-cell acute lymphoblastic leukaemia (B-ALL) and large deletions of 6q [del(6q)] as a model, we integrated analysis of functional cDNA clone tracking assays with patient genomic and transcriptomic data, to identify the transcription factors FOXO3 and PRDM1 as candidate tumour suppressor genes (TSG). Analysis of cell cycle and transcriptomic changes following overexpression of FOXO3 or PRDM1 indicated that they co-operate to promote cell cycle exit at the pre-B cell stage. FOXO1 abnormalities are absent in B-ALL, but like FOXO3, FOXO1 expression suppressed growth of TCF3::PBX1 and ETV6::RUNX1 B-ALL in-vitro. While both FOXOs induced PRDM1 and other genes contributing to late pre-B cell development, FOXO1 alone induced the key transcription factor, IRF4, and chemokine, CXCR4. CRISPR-Cas9 screening identified FOXO3 as a TSG, while FOXO1 emerged as essential for B-ALL growth. We relate this FOXO3-specific leukaemia-protective role to suppression of glycolysis based on integrated analysis of CRISPR-data and gene sets induced or suppressed by FOXO1 and FOXO3. Pan-FOXO agonist Selinexor induced the glycolysis inhibitor TXNIP and suppressed B-ALL growth at low dose (ID50 < 50 nM).

Whole genome sequencing provides comprehensive genetic testing in childhood B-cell acute lymphoblastic leukaemia.

Childhood B-cell acute lymphoblastic leukaemia (B-ALL) is characterised by recurrent genetic abnormalities that drive risk-directed treatment strategies. Using current techniques, accurate detection of such aberrations can be challenging, due to the rapidly expanding list of key genetic abnormalities. Whole genome sequencing (WGS) has the potential to improve genetic testing, but requires comprehensive validation. We performed WGS on 210 childhood B-ALL samples annotated with clinical and genetic data. We devised a molecular classification system to subtype these patients based on identification of key genetic changes in tumour-normal and tumour-only analyses. This approach detected 294 subtype-defining genetic abnormalities in 96% (202/210) patients. Novel genetic variants, including fusions involving genes in the MAP kinase pathway, were identified. WGS results were concordant with standard-of-care methods and whole transcriptome sequencing (WTS). We expanded the catalogue of genetic profiles that reliably classify PAX5alt and ETV6::RUNX1-like subtypes. Our novel bioinformatic pipeline improved detection of DUX4 rearrangements (DUX4-r): a good-risk B-ALL subtype with high survival rates. Overall, we have validated that WGS provides a standalone, reliable genetic test to detect all subtype-defining genetic abnormalities in B-ALL, accurately classifying patients for the risk-directed treatment stratification, while simultaneously performing as a research tool to identify novel disease biomarkers.

Integrative genomic analysis of childhood acute lymphoblastic leukaemia lacking a genetic biomarker in the UKALL2003 clinical trial.

Incorporating genetics into risk-stratification for treatment of childhood B-progenitor acute lymphoblastic leukaemia (B-ALL) has contributed significantly to improved survival. In about 30% B-ALL (B-other-ALL) without well-established chromosomal changes, new genetic subtypes have recently emerged, yet their true prognostic relevance largely remains unclear. We integrated next generation sequencing (NGS): whole genome sequencing (WGS) (n = 157) and bespoke targeted NGS (t-NGS) (n = 175) (overlap n = 36), with existing genetic annotation in a representative cohort of 351 B-other-ALL patients from the childhood ALL trail, UKALL2003. PAX5alt was most frequently observed (n = 91), whereas PAX5 P80R mutations (n = 11) defined a distinct PAX5 subtype. DUX4-r subtype (n = 80) was defined by DUX4 rearrangements and/or ERG deletions. These patients had a low relapse rate and excellent survival. ETV6::RUNX1-like subtype (n = 21) was characterised by multiple abnormalities of ETV6 and IKZF1, with no reported relapses or deaths, indicating their excellent prognosis in this trial. An inferior outcome for patients with ABL-class fusions (n = 25) was confirmed. Integration of NGS into genomic profiling of B-other-ALL within a single childhood ALL trial, UKALL2003, has shown the added clinical value of NGS-based approaches, through improved accuracy in detection and classification into the range of risk stratifying genetic subtypes, while validating their prognostic significance.

Genetic and genomic analysis of acute lymphoblastic leukemia in older adults reveals a distinct profile of abnormalities: analysis of 210 patients from the UKALL14 and UKALL60+ clinical trials.

Despite being predominantly a childhood disease, the incidence of acute lymphoblastic leukemia (ALL) has a second peak in adults aged 60 years and over. These older adults fare extremely poorly with existing treatment strategies and very few studies have undertaken a comprehensive genetic and genomic characterization to improve prognosis in this age group. We performed cytogenetic, single nucleotide polymorphism (SNP) array and next-generation sequencing (NGS) analyses on samples from 210 patients aged ≥60 years from the UKALL14 and UKALL60+ clinical trials. BCR-ABL1-positive disease was present in 26% (55/210) of patients, followed by low hypodiploidy/near triploidy in 13% (28/210). Cytogenetically cryptic rearrangements in CRLF2, ZNF384 and MEF2D were detected in 5%, 1% and <1% of patients, respectively. Copy number abnormalities were common and deletions in ALL driver genes were seen in 77% of cases. IKZF1 deletion was present in 51% (40/78) of samples tested and the IKZF1plus profile was identified in over a third (28/77) of cases of B-cell precursor ALL. The genetic good-risk abnormalities high hyperdiploidy (n=2), ETV6-RUNX1 (no cases) and ERG deletion (no cases) were exceptionally rare in this cohort. RAS pathway mutations were seen in 17% (4/23) of screened samples. KDM6A abnormalities, including biallelic deletions, were discovered in 5% (4/78) of SNP arrays and 9% (2/23) of NGS samples, and represent novel, potentially therapeutically actionable lesions using EZH2 inhibitors. Outcome remained poor with 5-year event-free and overall survival rates of 17% and 24%, respectively, across the cohort, indicating a need for novel therapeutic strategies.

PCR amplicons identify widespread copy number variation in human centromeric arrays and instability in cancer.

Centromeric α-satellite repeats represent ~6% of the human genome, but their length and repetitive nature make sequencing and analysis of those regions challenging. However, centromeres are essential for the stable propagation of chromosomes, so tools are urgently needed to monitor centromere copy number and how it influences chromosome transmission and genome stability. We developed and benchmarked droplet digital PCR (ddPCR) assays that measure copy number for five human centromeric arrays. We applied them to characterize natural variation in centromeric array size, analyzing normal tissue from 37 individuals from China and 39 individuals from the US and UK. Each chromosome-specific array varies in size up to 10-fold across individuals and up to 50-fold across chromosomes, indicating a unique complement of arrays in each individual. We also used the ddPCR assays to analyze centromere copy number in 76 matched tumor-normal samples across four cancer types, representing the most-comprehensive quantitative analysis of centromeric array stability in cancer to date. In contrast to stable transmission in cultured cells, centromeric arrays show gain and loss events in each of the cancer types, suggesting centromeric α-satellite DNA represents a new category of genome instability in cancer. Our methodology for measuring human centromeric-array copy number will advance research on centromeres and genome integrity in normal and disease states.

Dynamic clonal progression in xenografts of acute lymphoblastic leukemia with intrachromosomal amplification of chromosome 21.

Intrachromosomal amplification of chromosome 21 is a heterogeneous chromosomal rearrangement occurring in 2% of cases of childhood precursor B-cell acute lymphoblastic leukemia. These abnormalities are too complex to engineer faithfully in animal models and are unrepresented in leukemia cell lines. As a resource for future functional and preclinical studies, we have created xenografts from the leukemic blasts of patients with intrachromosomal amplification of chromosome 21 and characterized them by in-vivo and ex-vivo luminescent imaging, flow immunophenotyping, and histological and ultrastructural analyses of bone marrow and the central nervous system. Investigation of up to three generations of xenografts revealed phenotypic evolution, branching genomic architecture and, compared with other B-cell acute lymphoblastic leukemia genetic subtypes, greater clonal diversity of leukemia-initiating cells. In support of intrachromosomal amplification of chromosome 21 as a primary genetic abnormality, it was always retained through generations of xenografts, although we also observed the first example of structural evolution of this rearrangement. Clonal segregation in xenografts revealed convergent evolution of different secondary genomic abnormalities implicating several known tumor suppressor genes and a region, containing the B-cell adaptor, PIK3AP1, and nuclear receptor co-repressor, LCOR, in the progression of B-cell acute lymphoblastic leukemia. Tracking of mutations in patients and derived xenografts provided evidence for co-operation between abnormalities activating the RAS pathway in B-cell acute lymphoblastic leukemia and for their aggressive clonal expansion in the xeno-environment. Bi-allelic loss of the CDKN2A/B locus was recurrently maintained or emergent in xenografts and also strongly selected as RNA sequencing demonstrated a complete absence of reads for genes associated with the deletions.

Characterisation of the genomic landscape of CRLF2-rearranged acute lymphoblastic leukemia.

Deregulated expression of the type I cytokine receptor, CRLF2, is observed in 5-15% of precursor B-cell acute lymphoblastic leukaemia (B-ALL). We aimed to determine the clinical and genetic landscape of those with IGH-CRLF2 or P2RY8-CRLF2 (CRLF2-r) using multiple genomic approaches. Clinical and demographic features of CRLF2-r patients were characteristic of B-ALL. Patients with IGH-CRLF2 were older (14 y vs. 4 y, P < .001), while the incidence of CRLF2-r among Down syndrome patients was high (50/161, 31%). CRLF2-r co-occurred with primary chromosomal rearrangements but the majority (111/161, 69%) had B-other ALL. Copy number alteration (CNA) profiles were similar to B-other ALL, although CRLF2-r patients harbored higher frequencies of IKZF1 (60/138, 43% vs. 77/1351, 24%) and BTG1 deletions (20/138, 15% vs. 3/1351, 1%). There were significant differences in CNA profiles between IGH-CRLF2 and P2RY8-CRLF2 patients: IKZF1 (25/35, 71% vs. 36/108, 33%, P < .001), BTG1 (11/35, 31% vs. 10/108, 9%, P =.004), and ADD3 deletions (9/19, 47% vs. 5/38, 13%, P =.008). A novel gene fusion, USP9X-DDX3X, was discovered in 10/54 (19%) of patients. Pathway analysis of the mutational profile revealed novel involvement for focal adhesion. Although the functional relevance of many of these abnormalities are unknown, they likely activate additional pathways, which may represent novel therapeutic targets.

EBF1-PDGFRB fusion in pediatric B-cell precursor acute lymphoblastic leukemia (BCP-ALL): genetic profile and clinical implications.

The EBF1-PDGFRB gene fusion accounts for <1% of B-cell precursor acute lymphoblastic leukemia (ALL) cases and occurs within the Philadelphia-like ALL subtype. We report 15 EBF1-PDGFRB-positive patients from childhood ALL treatment trials (ALL 97/99, UKALL 2003, UKALL 2011) in the United Kingdom. The fusion arose from interstitial deletion of 5q33 (n = 11), balanced rearrangement (n = 2), or complex rearrangement (n = 2). There was a predominance of females (n = 11), median age of 12 years, and median white blood cell count of 48.8 × 10(9)/L. Among 12 patients who achieved complete remission on earlier trials (ALL 97/99 and UKALL 2003), 10 were positive for minimal residual disease (MRD) at the end of induction, and 7 relapsed 18 to 59 months after diagnosis. The majority (9 of 12) remained alive 6 to 9 years after diagnosis. There are reports of EBF1-PDGFRB-positive patients who are refractory to conventional chemotherapy who achieve complete response when treated with the tyrosine kinase inhibitor imatinib. These findings have prompted screening for EBF1-PDGFRB in patients entered onto the current UKALL 2011 trial for whom induction therapy failed, who did not achieve remission by day 29, or who remained MRD positive (>0.5%) at week 14. Two UKALL 2011 patients, positive for EBF1-PDGFRB, received imatinib; 1 died 6 months after a matched unrelated bone marrow transplant as a result of undefined encephalopathy, and the other remained in remission 10 months after diagnosis.

Combined MYC and P53 defects emerge at medulloblastoma relapse and define rapidly progressive, therapeutically targetable disease.

We undertook a comprehensive clinical and biological investigation of serial medulloblastoma biopsies obtained at diagnosis and relapse. Combined MYC family amplifications and P53 pathway defects commonly emerged at relapse, and all patients in this group died of rapidly progressive disease postrelapse. To study this interaction, we investigated a transgenic model of MYCN-driven medulloblastoma and found spontaneous development of Trp53 inactivating mutations. Abrogation of p53 function in this model produced aggressive tumors that mimicked characteristics of relapsed human tumors with combined P53-MYC dysfunction. Restoration of p53 activity and genetic and therapeutic suppression of MYCN all reduced tumor growth and prolonged survival. Our findings identify P53-MYC interactions at medulloblastoma relapse as biomarkers of clinically aggressive disease that may be targeted therapeutically.

DNA methylation profiling of medulloblastoma allows robust subclassification and improved outcome prediction using formalin-fixed biopsies.

Molecular subclassification is rapidly informing the clinical management of medulloblastoma. However, the disease remains associated with poor outcomes and therapy-associated late effects, and the majority of patients are not characterized by a validated prognostic biomarker. Here, we investigated the potential of epigenetic DNA methylation for disease subclassification, particularly in formalin-fixed biopsies, and to identify biomarkers for improved therapeutic individualization. Tumor DNA methylation profiles were assessed, alongside molecular and clinical disease features, in 230 patients primarily from the SIOP-UKCCSG PNET3 clinical trial. We demonstrate by cross-validation in frozen training and formalin-fixed test sets that medulloblastoma comprises four robust DNA methylation subgroups (termed WNT, SHH, G3 and G4), highly related to their transcriptomic counterparts, and which display distinct molecular, clinical and pathological disease characteristics. WNT patients displayed an expected favorable prognosis, while outcomes for SHH, G3 and G4 were equivalent in our cohort. MXI1 and IL8 methylation were identified as novel independent high-risk biomarkers in cross-validated survival models of non-WNT patients, and were validated using non-array methods. Incorporation of MXI1 and IL8 into current survival models significantly improved the assignment of disease risk; 46 % of patients could be classified as 'favorable risk' (>90 % survival) compared to 13 % using current models, while the high-risk group was reduced from 30 to 16 %. DNA methylation profiling enables the robust subclassification of four disease subgroups in frozen and routinely collected/archival formalin-fixed biopsy material, and the incorporation of DNA methylation biomarkers can significantly improve disease-risk stratification. These findings have important implications for future risk-adapted clinical disease management.

MYC family amplification and clinical risk-factors interact to predict an extremely poor prognosis in childhood medulloblastoma.

The MYC oncogenes are the most commonly amplified loci in medulloblastoma, and have previously been proposed as biomarkers of adverse disease prognosis by us and others. Here, we report focussed and comprehensive investigations of MYCC, MYCN and MYCL in an extensive medulloblastoma cohort (n = 292), aimed to define more precisely their biological significance and optimal clinical application to direct improved disease risk-stratification and individualisation of therapy. MYCC and MYCN expression elevations were multifactorial, associated with high-risk (gene amplification, large-cell/anaplastic pathology (LCA)) and favourable-risk (WNT/SHH molecular subgroups) disease features. Highly variable cellular gene amplification patterns underlay overall MYC copy number elevations observed in tumour biopsies; we used these alternative measures together to define quantitative methodologies and thresholds for amplification detection in routinely collected tumour material. MYCC and MYCN amplification, but not gain, each had independent prognostic significance in non-infants (≥3.0-16.0 years), but MYCC conferred a greater hazard to survival than MYCN when considered across this treatment group. MYCN's weaker group-wide survival relationship may be explained by its pleiotropic behaviour between clinical disease-risk groups; MYCN predicted poor prognosis in clinical high-risk (metastatic (M+) or LCA), but not standard-risk, patients. Extending these findings, survival decreased in proportion to the total number of independently significant high-risk features present (LCA, M+ or MYCC/MYCN amplification). This cumulative-risk model defines a patient group characterised by ≥2 independent risk-factors and an extremely poor prognosis (<15% survival), which can be identified straightforwardly using the reported MYC amplification detection methodologies alongside clinical assessments, enabling targeting for novel/intensified therapies in future clinical studies.

Rapid diagnosis of medulloblastoma molecular subgroups.

PURPOSE: Microarray studies indicate medulloblastoma comprises distinct molecular disease subgroups, which offer potential for improved clinical management. EXPERIMENTAL DESIGN: Minimal mRNA expression signatures diagnostic for the Wnt/Wingless (WNT) and Sonic Hedgehog (SHH) subgroups were developed, validated, and used to assign subgroup affiliation in 173 tumors from four independent cohorts, alongside a systematic investigation of subgroup clinical and molecular characteristics. RESULTS: WNT tumors [12% (21/173)] were diagnosed >5 years of age (peak, 10 years), displayed classic histology, CTNNB1 mutation (19/20), and associated chromosome 6 loss, and have previously been associated with favorable prognosis. SHH cases [24% (42/173)] predominated in infants (<3 years) and showed an age-dependent relationship to desmoplastic/nodular pathology; all infant desmoplastic/nodular cases (previously associated with a good outcome) were SHH-positive, but these relationships broke down in noninfants. PTCH1 mutations were common [34% (11/32)], but PTCH1 exon1c hypermethylation, chromosome 9q and REN (KCTD11) genetic loss were not SHH associated, and SMO or SUFU mutation, PTCH1 exon1a or SUFU hypermethylation did not play a role, indicating novel activating mechanisms in the majority of SHH cases. SHH tumors were associated with an absence of COL1A2 methylation. WNT/SHH-independent medulloblastomas [64% (110/173)] showed all histologies, peaked at 3 and 6 years, and were exclusively associated with chromosome 17p loss. CONCLUSIONS: Medulloblastoma subgroups are characterized by distinct genomic, epigenomic and clinicopathologic features, and clinical outcomes. Validated array-independent gene expression assays for the rapid assessment of subgroup affiliation in small biopsies provide a basis for their routine clinical application, in strategies including molecular disease-risk stratification and delivery of targeted therapeutics.

Definition of disease-risk stratification groups in childhood medulloblastoma using combined clinical, pathologic, and molecular variables.

PURPOSE: Medulloblastomas are heterogeneous and include relatively good-prognosis tumors characterized by Wnt pathway activation, as well as those that cannot be successfully treated with conventional therapy. Developing a practical therapeutic stratification that allows accurate identification of disease risk offers the potential to individualize adjuvant therapy and to minimize long-term adverse effects in a subgroup of survivors. METHODS: Using formalin-fixed paraffin-embedded (FFPE) tissue for immunohistochemistry, fluorescent in situ hybridization, and direct sequencing to identify tumors with a Wnt pathway signature and those harboring copy number abnormalities (CNAs) of potential prognostic significance (MYC/MYCN amplification, CNAs of chromosome 6 and 17), we evaluated clinical, pathologic, and molecular outcome indicators and stratification models in a cohort (n = 207) of patients with medulloblastoma 3 to 16 years of age from the International Society of Pediatric Oncology CNS9102 (PNET3) trial. RESULTS: Metastatic disease and large-cell/anaplastic (LC/A) phenotype were the clinicopathologic variables associated with poor progression-free survival (PFS). Nuclear immunoreactivity for ß-catenin, CTNNB1 mutation, and monosomy 6 all identified a group of good-prognosis patients. MYC amplification was associated with poor outcome, but other CNAs were not. Low-risk medulloblastomas were defined as ß-catenin nucleopositive tumors without metastasis at presentation, LC/A phenotype, or MYC amplification. High-risk medulloblastomas were defined as tumors with metastatic disease, LC/A phenotype, or MYC amplification. Low-risk, standard-risk, and high-risk categories of medulloblastoma had significantly (P < .0001) different outcomes. CONCLUSION: Integrating assays of molecular biomarkers undertaken on routinely collected diagnostic FFPE tissue into stratification schemes for medulloblastoma alongside clinical and pathologic outcome indicators can refine current definition of disease risk and guide adjuvant therapy.

Amplification and overexpression of Hsa-miR-30b, Hsa-miR-30d and KHDRBS3 at 8q24.22-q24.23 in medulloblastoma.

BACKGROUND: Medulloblastoma is the most common malignant brain tumour of childhood. The identification of critical genes involved in its pathogenesis will be central to advances in our understanding of its molecular basis, and the development of improved therapeutic approaches. METHODOLOGY/PRINCIPAL FINDINGS: We performed a SNP-array based genome-wide copy number analysis in medulloblastoma cell lines, to identify regions of genomic amplification and homozygous deletion, which may harbour critical disease genes. A series of novel and established medulloblastoma defects were detected (MYC amplification (n = 4), 17q21.31 high-level gain (n = 1); 9p21.1-p21.3 (n = 1) and 6q23.1 (n = 1) homozygous deletion). Most notably, a novel recurrent region of genomic amplification at 8q24.22-q24.23 was identified (n = 2), and selected for further investigation. Additional analysis by interphase fluorescence in situ hybridisation (iFISH), PCR-based mapping and SNP-array revealed this novel amplification at 8q24.22-q24.23 is independent of MYC amplification at 8q24.21, and is unique to medulloblastoma in over 800 cancer cell lines assessed from different tumour types, suggesting it contains key genes specifically involved in medulloblastoma development. Detailed mapping identified a 3Mb common minimal region of amplification harbouring 3 coding genes (ZFAT1, LOC286094, KHDRBS3) and two genes encoding micro-RNAs (hsa-miR-30b, hsa-miR-30d). Of these, only expression of hsa-miR-30b, hsa-miR-30d and KHDRBS3 correlated with copy number status, and all three of these transcripts also displayed evidence of elevated expression in sub-sets of primary medulloblastomas, measured relative to the normal cerebellum. CONCLUSIONS/SIGNIFICANCE: These data implicate hsa-miR-30b, hsa-miR-30d and KHDRBS3 as putative oncogenic target(s) of a novel recurrent medulloblastoma amplicon at 8q24.22-q24.23. Our findings suggest critical roles for these genes in medulloblastoma development, and further support the contribution of micro-RNA species to medulloblastoma pathogenesis.